Forces exerted by a correlated fluid on embedded inclusions 
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We investigate the forces exerted on embedded inclusions by a fluid medium with long-range 
correlations, described by an effective scalar field theory. Such forces are the basis for the medium- 
mediated Casimir-like force. To study these forces beyond thermal average, it is necessary to define 
them in each microstate of the medium. Two different definitions of these forces are currently used 
in the literature. We study the assumptions underlying them. We show that only the definition 
that uses the stress tensor of the medium gives the sought-after force exerted by the medium on an 
embedded inclusion. If a second inclusion is embedded in the medium, the thermal average of this 
force gives the usual Casimir-like force between the two inclusions. The other definition can be used 
in the different physical case of an object that interacts with the medium without being embedded 
in it. We show in a simple example that the two definitions yield different results for the variance 
of the Casimir-like force. 



I. INTRODUCTION AND SUMMARY 

The present work deals with the force exerted by a 
fluid medium with long-range correlations, described by 
an effective scalar field theory, on embedded inclusions. 
As inclusions constrain the fluctuations of the correlated 
fluid medium, the force exerted by the medium on one 
inclusion depends on the position of the other inclusions. 
Thus, medium-mediated interactions appear between in- 
clusions. These fluctuation-induced interactions are anal- 
ogous to the Casimir force, which arises between two 
uncharged metallic plates in a vacuum because of the 
boundary conditions imposed by the plates on the quan- 
tum fluctuations of the electromagnetic field [1]. Casimir- 
like effects driven by thermal fluctuations of fluid media 
were first discussed by Fisher and de Gennes in the con- 
text of critical mixtures [2] , and they have recently been 
measured [3]. Such forces appear in many other systems, 
including liquid crystals, fluid membranes, fluid inter- 
faces and superfluids (for a review, see Refs. [4, 5]). 

The Casimir-like force between two inclusions is usu- 
ally defined as —dF/dl, where F is the free energy of 
the fluid medium with two inclusions separated by a dis- 
tance H. However, this definition only provides a thermal 
average force. In order to study the fluctuations of the 
force, as well as out-of-equilibrium situations, it is neces- 
sary to define the force exerted on an inclusion by a fluid 
medium in each microstate of this medium. Two differ- 
ent deflnitions of this force are currently used: the first 
one uses the stress tensor of the medium [6-12], while 
the second one is based on differentiating the effective 
Hamiltonian with respect to the position of the inclusion 
while keeping constant the field that describes the state 
of the medium [13-16]. The aim of the present work is to 
clarify the difference between these two definitions and 
to determine their respective domains of validity. 

In this paper, we consider a correlated fluid medium 
described by a scalar field 0: each microstate of this ef- 
fective field theory corresponds to a given 4>. We study 



the force / exerted by this fluid medium on an embed- 
ded inclusion. This force / is defined as the negative 
gradient of the effective Hamiltonian H of the medium 
together with the inclusion, with respect to the position 
of the inclusion. The validity of this fundamental defini- 
tion in our coarse-grained description is justified in Ap- 
pendix A. In order to determine / in a given microstate 
of the medium, the gradient of H must be taken in this 
microstate. There are two distinct ways of interpreting 
these words, yielding two different routes to calculate the 
force. In the first route, the field is kept constant in 
the Eulerian sense, i.e., (j) remains the same at each point 
in space. In the second route, is kept constant in the 
Lagrangian sense, i.e., each fluid particle of the medium 
keeps the same value of cj) during the infinitesimal defor- 
mation that is associated with the displacement of the 
inclusion. We show that the second route gives the inte- 
gral of the stress tensor of the medium on the boundary 
of the inclusion, which corresponds to the definition used 
in Refs. [6-12]. We argue that this definition is the right 
one for an embedded inclusion. We also consider the case 
of non-embedded influencing objects, which interact with 
the medium without being embedded in it, and we argue 
that the first route, which corresponds to the definition 
used in Refs. [13-16], is the right one in such situations. 

This paper is organized as follows: in Sec. II, we study 
the force / exerted on an inclusion by the fluid medium 
in a given microstate, starting from the variation of the 
total energy of the system when it undergoes a generic 
infinitesimal deformation, and we compare the two routes 
introduced above. In Sec. HI, we proceed similarly in the 
case of a non-embedded influencing object. In Sec. IV, we 
show the link between the thermal average of the force / 
and the Casimir-like force. In Sec. V, we study a simple 
example of Casimir-like force, where we show that the 
variance of the force depends strongly on the route that 
is chosen. This example illustrates the importance of 
choosing the right definition when studying Casimir-like 
forces beyond their average value at thermal equilibrium. 
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Finally, Sec. VI is a conclusion. 



II. EMBEDDED INCLUSION 

Let us consider an infinite d-dimensional fluid medium 
(d > 1) with short-range interactions, described in a 
coarse-grained fashion by a scalar field cj) with Hamilto- 
nian density V0). Let us assume that an embedded 
inclusion with center of mass in £ G M.'^ extends over the 
region A of this medium (see Fig. 1). The inclusion is 
made of a different material, typically a solid, where the 
physical field (j) is not defined. For instance, in critical bi- 
nary mixtures, is the order parameter, i.e., the shifted 
concentration in one component: cj) = c—c"^"^, where c'^'''' 
is the critical concentration of this component [10], so (j) 
is not defined in solid objects immersed in the mixture. 




FIG. 1. (Color online) Schematic representation of an in- 
clusion (white) embedded in a two-dimensional fluid medium 
(shaded). The inclusion extends over A, which is centered on 
so (f) is defined only on B = \ ^. The contour of A is 
called 5, and n denotes its outward normal. 

Since the embedded inclusion extends over the region 
A, delimited by the hypersurface S, (j) is defined only in 
the region B ^ \ A (see Fig. 1). We assume that 
the field (j> is coupled to the inclusion at the interface S 
through short-range interactions modeled by a potential 
y ((/)), so the effective Hamiltonian of the system reads 



H = 



(1) 



where i?inci represents the internal energy of the inclu- 
sion, which will not be varied in the following. Note that 
our work can be generalized to a more general coupling 
potential, such as V{(j),^(j)). While our proofs will be 
presented with V{(j)) for the sake of simplicity, the final 
results for the force with 1^(0, V0) will be mentioned. 
Our aim is to calculate the force /, which is defined as 



dH 



(2) 



where is a unit vector in the i direction. The Einstein 
summation convention has been used in Eq. (2), and will 



be used throughout. The validity of the fundamental def- 
inition Eq. (2) in our coarse-grained description is justi- 
fied in Appendix A from the principle of virtual work. 

In order to obtain /, we must calculate the variation 
of H when the inclusion undergoes the infinitesimal dis- 
placement d£. Since it is impossible to move an embed- 
ded inclusion without moving the fiuid particles of the 
surrounding medium, we will consider a displacement 
field dr in the whole medium. We use the expression 
"fiuid particle" in a similar way as in fluid mechanics [17] : 
it denotes a (macroscopically) closed mesoscopic part of 
the fiuid medium, the state of which is described by the 
value of the coarse-grained field 0. Let us consider the 
generic infinitesimal transformation 



0(r) 



dr{r) , 
) + ^0(t-) 



(3) 



More explicitly, the particle initially in r undergoes the 
infinitesimal displacement Sr{r), and the value of at a 
fixed point r in space is modified by the quantity 5(f>{r). 
The functions 6r and 6(j), defined respectively on M.'^ and 
on B, are assumed to be regular and to take small val- 
ues, of order e (i.e., 6(j) — 0{e) and \Sr\/L — 0{e) where 
L is the characteristic length of interest, e.g., the dis- 
tance between two inclusions in a study of the Casimir- 
like force). During this transformation, the total (La- 
grangian) variation of 4> for the fluid particle initially in r 
is (5^0(r) = {(j>+S(j>){r+dr)-(l){r) = S(l){r)+V(j>{r)-dr{r) 
at first order in e. 

Let us assume that in the region A, Sr is constant, 
equal to di: thus, the inclusion undergoes an infinitesi- 
mal global translation that does not affect its internal en- 
ergy i?inci- The total variation of the energy then reads, 
at first order in e: 



6H 



dU ^ f dU 

Oi 



an . 



nse. 



(jxTr 



8V 



(4) 



where we have introduced n, the exterior normal to S 
(see Fig. 1), and we have used the notation 9^0 = d<p/dri 
for i € {1, . . . , d}. In this equation, the hypervolume in- 
tegral on B contains the standard Euler-Lagrange term. 
Meanwhile, the first term of the first hypersurface inte- 
gral comes from 



Hd^r 



(5) 



ISB Js 

while its second term is obtained via Stokes' theorem: 
dU _\ ,^ f dU 



d,. 



B 



d'W = - 



5(j)n,d'^-^r . (6) 



We may now calculate the force /, as defined in Eq. (2), 
in a given microstate of the fluid medium. We will ex- 
amine successively the two different routes presented in 
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the introduction, which correspond to two different ways 
of keeping (j) constant. 



A. First route 

The first way we may proceed is to keep <f) constant at 
each point in space during the infinitesimal transforma- 
tion, or, in other words, to keep the field (j) constant in 
the Eulerian sense. This amounts to taking 5^ = in 
Eq. (4), which gives: 



(1) 



dH 



,Eul. 



dV ^ 



r Bi 



(7) 

In the case where the coupling potential is V{(l), ^(j)) 
instead of a term —{dV / d{dj(j)))didj4> has to be 

added in the brackets in Eq. (7). 

However, it is physically not clear why each point in 
space should keep the same value of </) during a displace- 
ment in which each fluid particle of the system moves by 
5r. Another more formal argument also shows that it is 
artificial to keep the function (f) constant while moving 
infinitesimally the inclusion: the domain of definition of 
(j) itself dc;pc!nds on the position of the inclusion. For in- 
stance, when the inclusion is moved from A (centered on 

to A' (centered on + 51), the initial (j) is not defined 
in the region ^ \ where it should exist after the trans- 
formation. One way to deal with this mathematical issue 
is to consider the analytic continuation of (j) in this small 
region, and to truncate ^ in ^' \ A, but the physical 
meaning of this process is unclear. 

In other words, this first route is not adapted to cal- 
culate the force on an embedded inclusion because the 
position £ of the inclusion and the Eulerian field (p are 
not independent variables. Let us now move on to the 
second route. 



B. Second route 

Let us consider the example of critical binary mixtures, 
where (f) is the (shifted) concentration: during a displace- 
ment that is smooth at the microscopic scale, each fluid 
particle keeps the same concentration, so (j) is transported 
by fluid particles. Similarly, in the case of liquid crys- 
tals, the order parameter field is transported. The case 
of membranes and interfaces is a little bit more complex 
(see the discussion in Sec. II C). 

Let us focus on fields that are transported by fluid par- 
ticles during a deformation. For such a field 0, the correct 
force / will be provided by the second route defined in the 
introduction, where (j) is kept constant in the Lagrangian 
sense, i.e., each fiuid particle of the medium keeps a con- 
stant (j) during the infinitesimal transformation. In this 
case, 6'^<f) = 0, so that (5</)(r) = — V(/)(r) • 6r{r) for all r 



in B: Eq. (4) becomes 

51 



5H = - 



on 

5(5,<^) 
dH 



di65ri<fr 



Hi d' 



id-l„ 



(8) 



Using the stress tensor T of the fluid medium, which is 
discussed and derived in Appendix B, and which reads 



U6i 



d{dj4>) 



(9) 



and its divergence 



8 T 



9^ » 



dn 



dicf> = —di4>, (10) 



djTij Sn d' 



(11) 



we can rewrite Eq. (8) as 
5H = - [ 

JB JS 

This relation being valid for any continuous displace- 
ment field 6r such that Vr £ A, Sr = 6(., it yields the 
force /(^^ exerted by the medium on the inclusion, and 
the hypervolume density q of internal forces in the fiuid 
medium. Indeed, we can carry out an identification with 
Eq. (AlO), which comes from the principle of virtual work 
(see Appendix A). We obtain 



f(2) 



dH 



' d£i 

Vr e B, q{r) = 



, Lagr. 

5H 

Sri{r)' 



j TijHjd'^-^rei, (12) 



6j — ^j'^ijij^^ * 



(13) 



This second method of calculating / gives the integral 
of the stress tensor of the fluid medium on the boundary 
S of the inclusion. Contrary to f^^\ the force /"^^^ docs 
not depend on the short-range coupling V between the 
inclusion and the medium (in particular, Eq. (12) holds 
for V{4>, V(^) as well as for V{(j))). This comes from the 
continuity of the function Sr: the fluid particles of the 
nicidium that are infinitesimally close to the inclusion un- 
dergo the same small displacement as the inclusion, so the 
short-range interaction between the medium and the in- 
clusion does not vary during the displacement. In a real 
displacement in a viscous fluid, the displacement field 
has to be continuous at the inclusion boundary, and we 
consider virtual displacements consistent with this con- 
straint. Note that although the force /^^^ in each mi- 
crostate does not depend on V, the usual Casimir-like 
force, which is the thermal average of /^^^ (see Sec. IV), 
can depend on V since the statistical weight of each mi- 
crostate depends on V. 

Thus, the right deflnition of the force exerted by the 
medium on an embedded inclusion, obtained via the sec- 
ond route, is given by the integral of the stress tensor of 
the medium. This corresponds to the definition used in 
Refs. [6-12], and it is also very close to the definition used 
in Refs. [18, 19], which is based on integrating a density- 
dependent pressure on the boundary of the inclusion. 
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C. Domain of application 



D. Mean-field configuration 



We have considered a Hamiltonian density and a 
coupling potential V depending on (f> and V0. More gen- 
erally, H and V can depend on higher-order derivatives 
of (p. Our work can be adapted to such cases. In par- 
ticular, the force /^^^ exerted on an inclusion can be ex- 
pressed as the integral of the stress tensor of the medium 
also in these cases. The stress tensor associated with 
a Hamiltonian (or Lagrangian) density involving higher- 
order derivatives is discussed, e.g., in Refs. [20, 21]. 

We have focused on fields that are transported by fluid 
particles during a displacement, such as the order pa- 
rameter field in a critical mixture or in a liquid crystal. 
However, our reasoning can be adapted to other cases, 
especially when the field cj) has a direct relation to the 
position of the fluid particle. For instance, the Hamilto- 
nian density of a membrane described in the Monge gauge 
depends on the height <p{x, y) of the membrane with re- 
spect to a reference plane: for a displacement 5r{x, y) of 
the fluid particle initially in {x,y), the variation of the 
height (p at point {x,y) is 6(j) = Sr^ — dx(j)5rx — dyCpSvy 
instead of Sep = —di(j)Sri. A stress tensor can be defined 
for the membrane, taking this particularity into account 
[22-24]. Note that lipid membrane models (e.g., the Hel- 
frich model [25]) involve the membrane curvature, and 
thus, second derivatives of the height (p. We showed in 
Ref. [11] that the thermal average of the integral of the 
membrane stress tensor around a pointlike inclusion, in 
presence of a second inclusion, gives the usual Casimir- 
like force between the two inclusions, and we also used 
the membrane stress tensor to study the fluctuations of 
this force. 

In spite of this wider domain of application, the present 

work is restricted to the case of a fluid medium described 
by a scalar field (p such that its change during an in- 
flnitesimal displacement is a function of this infinitesimal 
displacement. More precisely, the relation S(p = —diCpSri 
(or its equivalent for the membrane) is crucial in our 
derivation of /'•^•'. As a counterexample, let us consider 
a one-dimensional elastic solid described by a scalar de- 
formation field (p{X) = x{X)—X, where X is the position 
of a particle in the reference (nondeformed) configuration 
of the solid, while x{X) is its position in the deformed 
configuration considered. In this case, changing the refer- 
ence coordinate X of a particle has no link with changing 
the deformation field (p. 

In addition, the assumption that the displacement 6r 
is continuous at the boundary of the inclusion is valid for 
a viscous fiuid, however small its viscosity, but not for an 
ideal fiuid. In the latter case, only the component of the 
displacement normal to the inclusion boundary has to be 
continuous. Thus, our approach may not be fully adapted 
to a superfluid. However, it is adapted to study Casimir- 
like forces in a superfiuid in the parallel plate geometry, 
which is the one usually considered. Indeed, modifying 
the distance between the plates involves a displacement 
perpendicular to the plates. 



The quantities f^^^ and /^^' obtained via the two dif- 
ferent routes are, in general, not equal. However, they 
cannot be distinguished in the most probable configura- 
tion of the field (p, i.e., in the mean-field configuration. 

Indeed, this configuration is such that the energy H is 
stationary with respect to variations of (p: 



= 5H 



d{d,(P 



+ 



L 



dV 



dn 



dcp d{di(P) 



(14) 



for any 54>. The bulk equilibrium condition yields the 
usual Euler-Lagrange equation, valid on B: 



dn 

d{di<p 



0, 



(15) 



while the boundary equilibrium condition gives the fol- 
lowing relation, valid on S: 



dV 



dn 

d{di<p 



-rii . 



(16) 



When Eq. (16) holds, /^-^^ and /^^^ are identical. Thus, 
the difference between the two routes is irrelevant when 
one considers the mean-field configuration. 



III. NON-EMBEDDED INFLUENCING OBJECT 

Let us now study the case of an object that interacts 
with the fluid medium without being embedded in it (see 
Fig. 2). For instance, it may be an optical trap creating 
a local electromagnetic field in the medium, or a protein 
binding very softly onto a lipid membrane so that the 
latter keeps its fiuidity: in these cases, there is no ma- 
terial object immersed in the medium. We will refer to 
this type of object as an "influencing object" , to make 
the distinction with the embedded inclusion. We will see 
that, contrary to the case of the embedded inclusion, it 
is here the first route that gives the correct force /. 

In the case of the influencing object, the fleld (p with 
Hamiltonian density n is defined everywhere in M** (re- 
call that the d-dimensional fluid medium is assumed to be 
inflnite). In the region A, the medium is affected by the 
influencing object: a term V{<p), representing the inter- 
action between the influencing object and the medium, 
adds to n. The effective Hamiltonian H of the fluid 
medium with the influencing object reads 



H 



I n{<p,v 

I [n + VlA\d'^r 

J 9."^ 



(fr+ I V{<p)d'^r 
A 



(17) 
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FIG. 2. (Color online) Schematic representation of a two- 
dimensional fluid medium with an influencing object. The 
object, represented by a black sphere, is not embedded in 
the medium (here, it is above the plane where the medium 
stands). It influences the medium, i.e., the field (j), in a region 
A, which is centered on t. The other notations are the same 
as in Fig. 1. 



where 1^ is the indicator function of the region A. Here 
too, our work can be generalized to 1^(0, V0), and the 
final results for the force will be mentioned in this case. 

In order to calculate the force / defined in Eq. (2), 
we can follow the same path as in the case of the embed- 
ded inclusion, by applying the generic transformation (3). 
The variation of H during this transformation reads, at 
first order: 



5H 



^S(j)d'^r+ [ V{^)n,Se,d'^-\, (18) 
0<P Js 



where the functional derivative of H with respect to is 



6H__dn dV 
6(j) d<f) d(j) 



Ia - d,. 



dn 

9(9,0) 



(19) 



A. First route 

If (f) is kept constant in the Eulerian sense during the 
infinitesimal displacement (i.e., 6(j) = 0), we obtain 



f(i) - _e- — 



dH 



, Eul. 



V{(t>)nd'^-^r . (20) 



The result for V0) is exactly similar. 

In fact, since the influencing object is not embedded in 
the medium, the force exerted on it by the medium in a 
given microstate can be calculated directly by moving the 
object with respect to the medium in a given configuration 
of 0. In other words, t and the Eulerian field (f) can 
be considered as independent variables. Then, we only 
need to take into account the variation of the interaction 
energy Ep{£) — Jj^V{(f>) d'^r between the medium and 
the object when the position £ of the object is modified: 



, Eul. 



- / ¥{(/>) nd 
Js 



This derivation is more physical than considering the full 
variation of H during the generic transformation (3). 
It shows that here, j'-'^' is simply the negative gradi- 
ent of the potential energy of interaction Ep between the 
medium and the influencing object. 

This definition of the force, which is the correct one for 
influencing objects that are not embedded in the medium, 
is the one used in Refs. [13-16]. In these works, this def- 
inition of the force is used to investigate Casimir-like in- 
teractions out of equilibrium [13, 14]. Note however that 
Casimir-like interactions are usually studied between em- 
bedded inclusions and not between non-embedded influ- 
encing objects. Besides, this definition has also been used 
to investigate drag forces in classical fields [15, 16]. In the 
latter works, the example studied in detail is a pointlike 
magnetic field moving through an Ising ferromagnet: it 
qualifies as an influencing object, since nothing material 
is embedded in the ferromagnet. However, the present 
work shows that the application of these results to the 
diffusion of inclusions embedded in membranes, which is 
discussed in Refs. [15, 161, should be taken with caution. 



Second route 



The Hamiltonian in Eq. (17) can be used to describe 
a medium coupled to an influencing object, but also a 
medium containing a "perturbative embedded inclusion" , 
i.e., an inclusion that is only slightly different from the 
surrounding medium. Indeed, in the latter case, is 
deflned everywhere in M**, and the energy density is per- 
turbed by an extra term V inside the inclusion. For in- 
stance, in lipid membranes, domains with a lipid com- 
position different from that of the rest of the membrane 
can be described as perturbative embedded inclusions, as 
well as similar structures in nematic liquid crystals. 

As for any other inclusion, it is meaningless to move 
a perturbative inclusion while keeping constant in the 
Eulerian sense, since moving the inclusion displaces the 
surrounding fluid particles. If we assume, as before, that 
(f) is transported by fluid particles, the second route is the 
right one to calculate the force exerted on a perturbative 
embedded inclusion. Let us now calculate this force. 

If the field is kept constant in the Lagrangian sense 
during the inflnitesimal transformation, i.e., S(f>(r) = 
— V(/)(r) • Sr{r) for all r in K'', Eq. (18) becomes 



6H 



SH 



8^4) 5n d'^r 



(21) ^{Ia^ ' Is ^^'^^ '^''"'^} ' ^^^^ 
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where, again, B = M"^ \ This yields, using Eq. (19), 



OH 



dt 
SH 



Lagr. 



[ V{cP) 
Js 



nd'^-h 



an dv_ / dn 



L 



V(l)d\ 



(23) 



Given the singularities in the energy density H + V1_a 
on S, the integral over A has to be calculated us- 
ing the proccdinc defined in Appendix A, Eq. (A15). 
Since {dV / d(j))'V (j) = VF, the two terms involving V 
in Eq. (23) cancel. Thus, using Eq. (10), we obtain 



f(2) 



/ djTijd'^rei 
J A 



Jjijnjd'' 



(24) 



where T is the stress tensor of the fluid medium. This 
result, which is independent of V, remains the same for 

mv0). 

The expression for /^^^ is the same here as for the em- 
bedded inclusion. In fact, a perturbative inclusion is a 
particular inclusion with £^inci = Xa(^ + ^) 'i'^''"^ which 
verifies SEind = during our transformation, and with- 
out any explicit boundary coiipling. 

Thus, while it is the first route that gives the correct 
force exerted by the medium on an influencing object, the 
second route is the right one in the case of a perturbative 
embedded inclusion. 



C. Mean-field configuration 

The two forces f^^^ and f^^^ are, in general, not equal, 
and we have seen that they are relevant to very diff'erent 
physical situations. However, as in the case of an em- 
bedded inclusion (see Sec. II D), these two quantities are 
equal in the mean-fleld configuration of the system. In- 
deed, in this configuration, 5H/5(j) = 0, so Eq. (23) shows 
that is identical to f^'^\ 



D. Formal relation between the two types of forces 

Independently of the physical interpretations of the 
forces /(I) and Eq. (23) gives a formal relation be- 
tween these two forces: 



-L 



1^ 



V(l)d\. 



(25) 



Given the singularities in the energy density H + V1^ 

on S, the integral over A has to be calculated using the 
procedure defined in Appendix A, Eq. (A15). 



Let us introduce the stress tensor T' of the compos- 
ite medium comprising the perturbative embedded in- 
clusion: the force density at each point of the medium, 
even inside the perturbative embedded inclusion, is q' = 
djT-jCi (see Appendix C). Thus, the basic definition 
Eq. (All) of the force / enables us to write: 



f(2) 



/ q'{r)d'^r. 
J A 



(26) 



Here again, the integral over A has to be calculated using 
the procedure defined in Eq. (A15). Besides, if C is a 
region such that ^ C C, we have 

/ q'{r)d'^r= [ ^Vcj^d'^r- [ VU)nd'^-^r, (27) 
Jc Jc 0(p Js 

where we have used the expression of q' in Eq. (C2). 
Thus, we can write 



(1) 



q'{r)d\- [ ^V(j>d'' 
Jc 0(p 



(28) 



which corresponds exactly to formula (18) in [14]. How- 
ever, this equation relates f^^^ and f^^^ only if C = A: 
the formal relation between f^^^ and f^^^ is given by 
Eq. (25). 



IV. CASIMIR-LIKE FORCE 
A. Embedded inclusions 

Let us consider a fiuid medium comprising two em- 
bedded inclusions with respective centers of mass at the 
origin of the frame and at point £, and let us introduce 
the unit vector u = tjl. The Casimir-like force exerted 
on the inclusion in I by the other one is usually defined 
through jc = -udF/dt, where i^(£) = -fcBTln2'(£) is 
the free energy of the system. We are going to show that 
fc corresponds to the average at thermal equilibrium of 
the force /(^' exerted on the inclusion in I by the medium 
containing the other inclusion. 

Let us assume that the line joining the centers of mass 
of the two inclusions is a symmetry axis of the system. 
This assumption is valid in the standard case of parallel 
plates, as well as for pointlike and spherical inclusions. 
Then, the effective Hamiltonian H of the system only 
depends on and on the distance I between the two 
inclusions. At thermal equilibrium, the statistical weight 
of a configuration ([</>], ^) is given by e~^^^^'^^'^\ where 
the notation [4>\ represents a functional dependence on 
(j). At a given the partition function of the system is 



V(j) e 



(29) 



where the functional integral runs over the functions (j) 
defined in the domain S = R'' \ (^o U ), where At, is 
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the region where the inclusion with center of mass in t 
stands (and similarly for ^o)- 

Let us now introduce a replica of the system presented 
above, identical to it except that the center of mass of 
the second inclusion is at point ^ + The partition 

function of this replica reads 

Z{t + Se) = j V'4> e-/3«(['^l'^+«) . (30) 

Here, the functional integral runs over the functions (f> 
defined in the domain = M'' \ (AoUA^+seu)- This dif- 
ference with respect to Eq. (29) is symbolized by a prime 
on the functional measure. Let us consider a smooth in- 
vertible function r' : B ^ B', r i-^ r'{r) = r + dr{r), 
where the infinitesimal virtual displacement field 5r is 
such that = on the first inclusion while 6r = S£u 
on the second one, and |^r|/i? is small, say of order e. Let 
us assume that the function r' maps a field to a field 
^ = (j) + S<j) such that ^(r'(r)) = (j){r) for all r G B, and 
5(f) = —Sr-'V(p at first order in e. This is motivated by our 
assumption that each fluid particle of the medium keeps 
the same value of </> during any real smooth infinitesimal 
displacement. As this process can be inverted, r' yields a 
one-to-one mapping of the states ([(/)], ^) of the first sys- 
tem to the states {[(j) + d(j)],£ + 6£) of the replica. We may 
thus write 

V<t>=]ld[<t>{r)]= Yl d[4>{r'{r))]=V'4,, (31) 

reB r'{r)eB' 

where the continuous products must be understood as 

UreB'^i^l^i'^)] = limjv^-oo n^Li d[(t){rn)], where {r„} is a 
regular mesh of B [26] . Hence, Eq. (30) can be rewritten 
as 

Z{e + 5£) = j Vct> e-/3«([^+5^].^+«) , (32) 

so the difference of free energy between the replica and 
the original system reads at first order in e: 

5F = -knT —= I Vct> 5H = {6H) , (33) 

where (.) denotes the average at thermal equilibrium, 
while 5Z = Z{£ + S£) - Z{£), and 6H = H{[(j) + 6(f>],£ + 
S£)-Hm,£). 

The expression of SH is given by Eq. (11), so we obtain 

5F = - I (ft) Sn dfr -5£ [ {T^i)ni d'^-'r , (34) 
Jb Js 

where the axes have been chosen so that one of them is 
along u: TuiUi d'^~^r denotes the component along this 
axis of the force transmitted through the infinitesimal hy- 
persurface d'^~^r. Recall that q represents the hypcrvol- 
ume density of internal in the medium, given by Eq. (13). 

As no average external forces are imposed to the sys- 
tem, the hypervolume density of external forces w veri- 
fies {w) = 0. Neglecting inertia, which is possible if the 



timescales considered are sufficiently large, Newton's sec- 
ond law applied to each Qxiid particle of the medium gives 
q — —w (see Appendix A). Thus, we obtain {q) = 0. 
This relation can also be obtained formally: Eqs. (10) 
and (13) show that q = SH/Scj), and the thermal 
average of the latter quantity can be proved to vanish 
using the Schwinger-Dyson equation [14]. Thus, Eq. (34) 
becomes 

SF = -5£jj,Tui)ni d''" V = -5£ (Z^^)) • u , (35) 

where we have used the expression of in Eq. (12). 
Note that Eq. (35) is independent of the virtual displace- 
ment field 8r chosen at the beginning of this discussion to 
map the states of the system onto the state of its replica. 
We obtain from Eq. (35): 

fc = -^u=({f^^^)-u)u. (36) 

Since (/'^^•') is along u by symmetry, we can conclude 
that fc = 

B. Non-embedded influencing objects 

Although Casimir-like forces arc usually studied be- 
tween embedded inclusions, let us now consider non- 
embedded influencing objects. In this case, the partition 
function can also be written as 

Z{£) = jvcj, e-'^^d^l'^) , (37) 

but here the functional integral runs over the functions 
4> defined on W^. In contrast to the case of the embedded 
inclusions, £ can be varied at constant Eulerian field (j). 
Thus we can write directly 

d£ J ^ Z d£ \d£ / ' 

(38) 

Thus, in the case of influencing objects, it is the average 
of /(I) that gives -udF/d£. 

Some effective Hamiltonians, such as the one in 
Eq. (17) and the one studied in the next section, can de- 
scribe both influencing objects and embedded inclusions. 
For such effective Hamiltonians, our results show that 
the thermal average of / is the same in these two physi- 
cal cases, each of them being treated with the appropri- 
ate route. If there are two inclusions, this thermal aver- 
age force corresponds to the Casimir-like force —u dF/d£. 
Thus, the usual (i.e., thermal average) Casimir-like force 
does not depend on the route that is chosen. In spite 
of this degeneracy concerning thermal average, it is very 
important to distinguish the two physical cases as soon 
as one wishes to go beyond the thermal average force. 
This point will be illustrated in the next section. 
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V. A SIMPLE EXAMPLE OF CASIMIR-LIKE 
FORCE 

A. The situation 

Let us now present a simple example with a Hamil- 
tonian that can describe both influencing objects and 
embedded inclusions. In this example, we calculate the 
variance of the Casimir-like force, both in the usual case 
of embedded inclusions (via the second route) and in the 
particular case of non-embedded influencing objects (via 
the first route), and we find two very different results. 
Thus, this example illustrates the importance of the dis- 
tinction between the two physical cases. 

Let us consider an infinite one-dimensional fluid 
medium described by a dimensionless scalar field (j) with 
Hamiltonian density 



(39) 



where (/>' = d(j)/dx. The length scale L = y/n/m which 
appears in Ji corresponds to the correlation length of (6, 
as shown in the following. The energy Hq of the medium 
is such that 

pHo = l3 [ dx H{x) = \ I dxdy (l){x)0{x,y)ct>{y) , 

(40) 

where /3 = {k-B.T)~^ , and the symmetric operator O is 
defined by 



(3m - Pk--^ 
ax'' 



d{x - y) . 



(41) 



Let us assume that there arc two pointlikc inclusions or 
non-embedded influencing objects, in x = and x = £, 
where ^ > 0, and that their coupling to the field ^ is 
given by 



(42) 



Since is dimensionless, a is dimensionless too. When 
a — >■ 00, this quadratic coupling yields Dirichlet bound- 
ary conditions on the inclusions. Indeed, in this limit, 
the statistical weight e~'^(^o+^^) gf any configuration such 
that 7^ on an inclusion goes to zero. 

For the calculations that follow, let us assume that 
there is an external field u conjugate to </>: 



dxu{x)(j){x) . 



(43) 



In order to calculate the partition function of the system 
Z[u] = j'^(t> 6-^^"°+^+"'-'^ , (44) 



let us carry out a Hubbard-Stratonovich transformation 
[27]. Using the relation 



27ra 



d V exp 



R2 



-2^+i(t;i ct){0)+V2<f>{i)) 



(45) 

where {vi,V2) is a two-dimensional vector, and perform- 
ing the Gaussian integration on (p, we obtain 



Z\u\ = / d^v exp 



V 

2a 



dxdy S{x)Gix,y)Siy) 



(46) 



In this expression, Zq is the partition function of the 
medium with no inclusion or other influencing object, 
while S{x) = u{x) + iviS{x) + iv2S{x — (), and G is 
the Green function of O. The latter can be obtained 
using a Fourier transform since the medium is infinite 
and translation-invariant: 



G{x,y) = 



far 

knTL 
2k 



dq ■ 



Jq{x-y) 



exp - 



m 



\x-y\ 



(47) 



where L = \J K/m. In the absence of inclusions or other 
influencing objects, the correlation function of </> is G, so 
L represents the correlation length of 0. Therefore, we 
expect Casimir-like forces to be most important when the 
distance i between the objects is such that £ <C i. 

Performing the Gaussian integration on {vi,V2) in 
Eq. (46) yields 



Z[u] 



exp 



(2 



dx dy u{x) C{x, y) u{y) 



aVdet M 

(48) 

where C{x,y) = G{x,y) - A'^{x) M-^A{y), with 
A^ix) = {G{x,0) , G{x,e)) and 



/G(0,0) + a-i G{0,e) 
^^^-y G{0,1) Gii,i) + a-^ 



(49) 



It is straightforward to deduce the moments of the Gaus- 
sian variable (j){x) from Eq. (48): 



{ct>ix)c(>{y)) = - 



1 SZ 
Z Su{x) 
1 d^Z 
Z 6u{x) 5u{y) 



0]=0, (50) 

[u = 0\ = C{x,y). 

(51) 



Thus, the correlation function of the fleld (j) in the pres- 
ence of the two objects is 

{m<t>{y))-{m){m) = c{x,y) 

= G{x,y)-A^{x)M-^A{y). 

(52) 
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B. Average Casimir-like force 

The average Casimir-like force (/) between the two 
inclusions can be calculated by differentiating the free 
energy F = —ksT \n{Z[u = 0]) of the system, which can 
be obtained from Eq. (48): 



dF _ kBT djlndetM) 
"di ~ 2 d£ 



(fcBT)2a2i2 (^i_e-¥) +4kBTaKL + AK'^ 



(53) 



Taking the hard-constraint limit a — > oo, which will 
be denoted in the following by a subscript index h, we 
obtain: 



{h} = 



kBT 



(l-e«) 



e<^L 



kBT 
' 21 



(54) 



As expected, the Casimir-like force vanishes when i:^ L: 
when the distance between the inclusions is much larger 
than the correlation length of cj), one inclusion cannot feel 
the effect of the other one. We have emphasized the "crit- 
ical regime" £ <^ L, where the Casimir-like force is most 
important: it has a simple power-law dependence. 

Meanwhile, in the perturbative limit a — >■ 0, which 
will be denoted in the following by a subscript index p, 
we obtain, to lowest order in a: 



{kBTfa^Le 
1^ 



{kBTfa^ (L - 2£) 



Up) = - 

(55) 

Here too, we have emphasized the critical regime £ <^ L: 
this time, the leading term is independent of £, so the 
force is nearly constant in this regime. 



C. Variance of the Casimir-like force 

We are now going to illustrate the difference between 
embedded inclusions and non-embedded influencing ob- 
jects by studying the variance of the Casimir-like force. 
The two routes to calculate the force / are detailed in 
the particular case of pointlike embedded inclusions or 
non-embedded influencing objects in Appendix D. 



1. Preliminary remark 

The variance we are going to study is that of the force 
/ defined at the coarse-graining level where the system 
is described by ([</>], ^). This force is averaged over the 
fundamental microstates "/x — )• {[4>],£)" that yield the 
coarse-grained field (j) and the inclusion position £. In- 
deed, it is shown in Appendix A that the force exerted 



by the medium on the inclusion in the fundamental mi- 

crostate /x verifies 



{fn)ij,- 



/([</>], ^ 



(56) 



where represents the statistical average at 

equilibrium over the fundamental microstates fi — t- 
{[(/)],£.). Therefore, the variance A^f of the force / is 
smaller than A^/^. 

The difference between A^/^ and A^f is due to fluc- 
tuation modes with wavelengths smaller than the cutoff 
a of the field </>: such modes are averaged in the coarse- 
graining procedure leading to /. Thus, we expect A^f to 
decrease if a increases. The small-wavelength fluctuation 
modes are generically the fastest to equilibrate, so above 
a certain time resolution, it is right to consider A^ f in- 
stead of A^/^, and for even longer time resolutions, it 
is appropriate to consider larger values of a. This sub- 
telty, which is linked to the coarse-grained nature of the 
force /, arises in all studies of fluctuations of Casimir-like 
forces [6, 11]. 



2. Embedded inclusions 

In the case of pointlike embedded inclusions, the force 
exerted on the inclusion in £ by the medium containing 

the other inclusion is obtained via the second route. It 
is given by Eq. (D6), adapted to d = 1: /f^) = T{£+) - 
T{£-). 

Let us first show explicitly in the present example that 
the thermal average of Z^^) gives back the Casimir force 
(/). Here, the stress tensor of the medium is 



dn 



m 



(57) 



and Eq. (51) enables to express its average as 

{T{X)) : 



^C,,(a;,a:) + yC(,T,.T), (58) 



where we have introduced the notation Cxy{x,y) = 
{dxdyC]{x,y). Such notations will be used in the follow- 
ing. This average can be calculated thanks to Eq. (47) 
and (52). Note that G^yix, x) has to be regularized using 
the short-distance cutoff a of the theory: 




+ m 
Arctan (^)' 



(59) 



We then recover the result obtained in Eq. (53): (/'^^) = 
{T{£+)) - {T{£-)) = (/). Note that the cutoff-dependent 
term Eq. (59) vanishes when one calculates the difference 
between the average stress on each side of the inclusion. 
This is reminiscent of the calculation of the (average) 
quantum Casimir force, either from the energy [1] or from 
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the radiation pressure [28], where two infinite quantities 
subtract to give a finite force. However, in the variance, 
this cutoff-dependent term will no longer vanish. 

We may now proceed to calculate the variance A^/^^^ 
of/(2): 

a2/(2) = K{r,r) + k{£+,£+) -2K{r,e+) , (60) 

where K{x,y) = {T{x)T{y)) - {T{x)) {T{y)). Given the 
expression of the stress tensor in Eq. (57), its correlation 
function K involves quartic terms in (p (or (/>'). Since 
4>{x) is a centered Gaussian variable, we may use Wick's 
theorem to calculate K. It yields: 



-Cli£-J+)-C^y {£-,£+)] 



(61) 



This variance can be calculated from Eqs. (47), (52) and 
(59). Note that Gxix,x) and Gy{x,x) can be obtained 
from a regularized Fourier transform: 



Ga;{x,x) = -Gy{x,x) 



kBT 
2w 



l/a 



dq 



iq 



i/a Kq^+m 



0. 



(62) 



The full analytical expression of A^/^^^ is quite heavy, 
so we will only present its limiting behaviors for large 
and small a, and wc will focus on the regime where a < 
£ <^ L, since the Casimir effect is strongest in the critical 
regime i <^ L. More precisely, in each limit, we carry 
out expansions in the small parameter a/L after setting 
£ = Ca where C is a constant. 

Taking the hard-constraint limit a — > oo, and then 
keeping only the leading order in a/L, we obtain: 



AVf =AV(^)t + AVf 



where the first term is independent of i and a: 



A2/(2)t 



while the second one arises from the inclusions: 



(63) 



(64) 



(65) 



If we consider inclusions imposing Dirichlet boundary 
conditions (i.e., (j){0) = (j){£) = 0) instead of inclusions 
imposing a potential V, we recover the result in Eqs. (63, 
64, 65) for the force variance. In this case too, the av- 
erage force can be obtained either directly from F or by 
using the stress tensor, and it gives back (fh) in Eq. (54). 

In the perturbative limit a — > 0, let us keep the two 
lowest orders in a: the leading term corresponds to the 
situation where there is no inclusion, a = 0, so the effect 
of the inclusions appears in the subleading term. We 
obtain 



2f(2)t 



AVft 



(66) 



where the first term is the one in Eq. (64) , still at leading 
order in a/L, while 



TT Ka 



-1 



(67) 



where wc have kept the two leading orders in a/L to 
see the ^-dependence. Note that these results can be 
obtained directly from a perturbative expansion of the 
stress tensor correlation function K at first order in (3V. 
Besides, using this perturbative method, it is straight- 
forward to study the more general problem of inclusions 
with a finite size s: wc have checked that it gives back 
the average force Eq. (55) and the variance Eqs. (66, 67) 
in the small-size limit. 

Both in the hard-constraint limit and in the pertur- 
bative limit (and also in the intermediary regime), the 
variance of the force features a term A2/(2)t^ expressed 
in Eq. (64), which is independent of i and a. It cor- 
responds to the variance of the zero-average fluctuating 
force exerted on a point of the medium, or on a single 
inclusion, by the rest of the medium, in the absence of 
any other inclusion. This term depends on the cutoff a, 
and it decreases when a increases, in agreement with the 
preliminary remark above. The other terms depend on £, 
and this dependence has the same origin as the Casimir- 
like interaction itself: it comes from the constraints im- 
posed by the inclusions on the fluctuations of (f>. They 
also depend on the intensity of the coupling constant a, 
because they come from the constraints imposed by the 
inclusions. The presence of a term independent of £ and 
of smaller ones that depend on £ is similar to the results 
of Refs. [6, 11]. 



3. Other influencing objects 

In the case of pointlikc influencing objects that arc not 
embedded in the fluid medium, the force exerted on the 
object in £ by the medium is obtained via the first route. 
It is given by Eq. (D3), for d = 1: 



f'^ = -ci>'{i)—m)) 



-akBT<j)'{£)<j){l). (68) 



Eq. (51) enables to express its average as 

l^f^^^) =-akBTCx{£,£), (69) 

which can be calculated from Eq. (47), (52) and (62), 
yielding once more the result obtained in Eq. (53): 
(./'^^■') = (./)• We have thus verified on the present exam- 
ple that all the definitions (i.e., the direct derivation from 
the free energy, and the two routes defined in the intro- 
duction) give the same result for the thermal average of 
the force. 

In order to calculate the variance of f^^\ we use Wick's 
theorem as above. It gives 



AV^') =a2(fcBr)' [G{£,£)G.y{£,£) + 2Cl{£,£)\ . 



(70) 
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This variance can be calculated from Eqs. (47), (52), (59) 
and (62). Here again, we will only present the limiting 
behaviors of A'^f^^\ in the same regimes as for the em- 
bedded inclusion. 

In the hard-constraint limit a — >■ oo, keeping the lead- 
ing order in a, and then keeping the leading order in a/L, 
we obtain: 



a2 .(1) _ {kBT)^a 



a 



(71) 



In the perturbative limit a ^ 0, if we keep the two 

lowest orders in a (the lowest order term docs not depend 
on so wc also keep the next one), we obtain: 



L £ n 
■- + - + - 

a a 4 



(72) 

where we have kept only the leading orders in a/L, and 
also the subleading one in the term to see the i- 
dependence. This result can also be found directly from 
a perturbative expansion of ((/^^^)^) at second order in 
pV. This variance diverges when L — > oo, but the corre- 
lation function itself diverges in this limit (see Eq. (47)). 

The present results are very different from the ones 
concerning the embedded inclusions. First, there is no 
term independent of the coupling constant a here, since, 
as explained above, the first route amounts to taking the 
derivative of the potential energy of interaction between 
the medium and the influencing object. Besides, here, 
the variance diverges as a in the hard-constraint limit 
a ^ oo, while it converged in the case of the embed- 
ded inclusion. Here follows a physical explanation of 
this divergence: for each microstate (f>, let r and s be 
such that 1^(^)1 = and |^(0)| = a". For a ^ oo, 
only the states such that r < —1/2 and s < —1/2 can 
keep a significant statistical weight: they are such that 
exp(— /3F) > 1/e, while for the others, exp(— /3y) — > 
exponentially when a — >■ oo. Let us consider the statisti- 
cally significant states with highest r, i.e., the ones such 
that r = —1/2: they will yield the largest forces f^^\ 
For these states, (f'-^'^)'^ ~ a {kBTf(t>''^{l). Since there is 
no particular constraint on (j)'{t), (/'•^•')^ typically scales 
as a in these states, which explains the behavior of the 
variance. Much more qualitatively, if a is infinite, the 
only allowed states are such that </>(£) = 0, but (/>'(^) is 
generically nonzero, so </>(£ -|- St) ^ 0: if it moved to this 
position, the object would have an infinite energy. Thus, 
it feels an infinite restoring force. While such contribu- 
tions can cancel in the average, they add in the variance, 
which thus diverges with a. 

In this simple example, we have found very different re- 
sults for the variance of the Casiniir-like force for embed- 
ded inclusions and for non-embedded influencing objects. 
It illustrates the importance of distinguishing these two 
physical cases and the associated routes for calculating 
the force. 

The out-of-equilibrium behavior of Casimir-like forces 

also depends strongly on the route that is taken. The ap- 
proach to equilibrium after a quench at = has been 



studied for plates imposing boundary conditions to the 
field 0, using the second route in Ref. [10], and using the 
first route in Ref. [14]. The two routes give very differ- 
ent results for the case where one plate imposes Dirichlet 
boundary conditions while the other one imposes Neu- 
mann boundary conditions. Indeed, in Ref. [10], the force 
exerted on one plate is not equal to the one exerted on 
the other plate during the relaxation (but both converge 
to the same equilibrium value), while in Ref. [14], these 
two forces are always equal. Note that obtaining different 
forces on each plate is not ruled out by the action-reaction 
principle, since the forces at stake are forces exerted by 
the medium containing one plate on the other plate, and 
not forces exerted by one plate on the other plate. 



VI. CONCLUSION 

We have investigated the force exerted on an inclusion 
by a fluid medium with long-range correlations, described 
by a scalar field in a coarse-grained theory. If a sec- 
ond inclusion is embedded in the medium, the thermal 
average of this force gives the Casimir-like force between 
the two inclusions. In order to go beyond the thermal 
average force, it is necessary to define properly the force 
/ exerted on an inclusion by the medium in each mi- 
crostate. In practice, one must take the negative gradient 
of the effective Hamiltonian with respect to the position 
of the inclusion in a given microstate of the medium. We 
have shed light onto the subtlety of this task, showing two 
routes that can be taken to calculate this gradient. In the 
first route, </) is kept constant in the Eulerian sense, while 
in the second one, (j) is kept constant in the Lagrangian 
sense. 

In the usual case of an embedded inclusion, the po- 
sition of the inclusion and the Eulerian field (p are not 
independent variables, so one should not take the first 
route. Indeed, when an inclusion is displaced infinitesi- 
mally, the surrounding fluid particles are displaced too. 
In many physical cases, (j) is transported by the fluid par- 
ticles of the medium during a displacement, so the second 
route is the correct one. It gives the integral of the stress 
tensor of the medium on the boundary of the inclusion. 

We have also considered the case of infiuencing objects 
that interact with the medium without being embedded 
in it. Contrary to inclusions, such objects can be moved 
with respect to the medium at a given Eulerian field </>. 
Thus, the first route is adapted to calculate the force 
exerted by the medium on an influencing object. 

In a nutshell, the two formal routes for calculating the 
force; / apply to two different physical cases. We have 
discussed the practical importance of this distinction. 
First, in the mean-fleld configuration, the two routes give 
the same result. Then, for effective Hamiltonians that 
can describe both physical cases, the two routes give the 
same thermal average of the force. However, the differ- 
ence between these two routes becomes crucial as soon 
as one wants to study this force beyond its mean-field 
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value and its thermal average. We have shown in a sim- 
ple example that they yield very different results for the 
variance of the Casimir-like force. Besides, comparing 
Refs. [10] and [14] shows that the out-of-equilibrium be- 
havior of Casimir-like forces also depends on the route 
that is taken. 

Two definitions of the force, corresponding to our two 
routes, are currently used to study the fluctuations of the 
Casimir-like force and its out-of-equilibrium behavior. 
Our work shows that the second route, which gives the 
integral of the stress tensor, should be used when study- 
ing Casimir-like forces between embedded inclusions, in 
agreement with Refs. [6-12]. In contrast, the first route, 
which is used in Refs. [13-16], should be reserved to the 
case of non-embedded influencing objects. 



Appendix A: Definition of the force / from the 
principle of virtual work 

The aim of the present Appendix is to provide a justifi- 
cation of the fundamental definition Eq. (2) in our coarse- 
grained, effective description. Let us consider, as in the 
main text, an infinite d-dimensional fluid medium with 
short-range interactions, and let us assume that an em- 
bedded inclusion extends over the region AcR'^ of this 
medium. Let us consider a "fundamental" microstate /i 
of the system constituted by tlie fluid medium with the 
inclusion, i.e., a microstate deflned not by the coarse- 
grained field (j) and by the position £ of the center of mass 
of the inclusion, but by all the underlying microscopic de- 
grees of freedom. Let us call the hypervolume density 
of forces exerted by the exterior on the system, and 
the hypervolume density of forces exerted on a piece of 
the system by the rest of the system, in the microstate /u. 
Both Wf_i and can take different values in the differ- 
ent fundamental microstates "/i — )■ ([0], €)" that yield the 
same coarse-grained field 4> and inclusion position £. The 
coarse-graining procedure regroups these microstates so 
that 



(Al) 



where E^j, is the energy of the fundamental microstate 
/U, while H{[(j)],t) is the effective energy of the coarse- 
grained microstate ([0],£). Let us now introduce the 
average q( [</>],€) (respectively, w{[(f)\,€)) of (respec- 
tively, Wij) over the microstates /i — >■ ([(/)],£): q and w 
are coarse-grained force densities. Explicitly, we have 



9p 



(A2) 



and similarly for w. 

Let us consider a material particle of hypervolume d'^r 
of the system: it can be a fluid particle of the medium 

in the sense defined in the main text, or a piece of the 
inclusion. Since it is a (macroscopically) closed particle. 



Newton's second law applies to it. In the fundamental 
microstate /x, it reads w^df^r + q^j.d'^r = d{pfj,d'^r)/dt, 
where is the hypervolume density of momentum in 
the microstate ji. This is equivalent to writing that ■ 
Sr d'^r + qfj, - Sr d'^r = [d(p^ d'^r) /dt] ■ 6r for any smooth 
infinitesimal virtual displacement field Sr. The former 
relation can be integrated on the whole system: 

/ w^-Srd''r+ [ q^-5rd'^r=j ^^El^^ ■ 5r . 

Jm'' jRti jR<i dt 

(A3) 

The first integral in Eq. (A4) represents the work SW^ of 
the external forces on the whole system, and is therefore 
equal to the variation SEf^ of the energy of the system 
during the infinitesimal displacement. Thus, averaging 
over the microstates /x ^ {[(f)], £), we obtain 



q{[<j)],£) ■drd'^r 



d{p^ d'^r) 



dt 



■ Sr , (A4) 



where we have used Eq. (A2). Neglecting the average 
inertial term in Eq. (A4) yields 



,£)-Srd\. 



(A5) 



It is possible to neglect the inertial term even out of equi- 
librium provided that the timescales considered are suf- 
ficiently large, as in case of the Langevin equation [29]. 

Let us consider the variation SH of the coarse-grained 
Hamiltonian H during the virtual displacement Sr. Let 
us denote hy 4> = (f> + S(f) the field after the displacement, 
and hy £ + S(. the position of the center of mass of the 
inclusion after the displacement. We can write at first 
order 



5H{[cP],£) 



H{[(l)],£ + S£)-H 

kBT ( ^ 



Z{[ct>],£) 



,£ + S£) 



(A6) 



We assume that the smooth infinitesimal virtual displace- 
ment Sr yields a one-to-one mapping of each fundamen- 
tal microstate ^ —>■([(/)],£) to a fundamental microstate 
/i' — )• ([^],€ -I- Si). In the present Appendix, we do not 
discuss the way cp is affected by the infinitesimal displace- 
ment — this problem is dealt with in the body of our 
paper — but we only assume that (f) is modified by an 
amount 5(f), which is a function of Sr. Then, 

Z{[^],£ + S£)= ^ e-^"^^' 



E 



-l3{E^+SE^) 



and we obtain at first order 



6Hm,£) = {SE^)^^^[^],e) 



(A7) 



(A8) 
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so Eq. (A5) becomes 

6Hm,i) = - [ qm,£)-Srd''r. (A9) 

Assuming tliat 6r is constant in A, equal to Si (so that 
the inclusion undergoes a translation), Eq. (A9) becomes 

SHm,i)^-f{[<P],e)-de- [ q{[<P],e)-Srd''r, (AlO) 

Jb 



where 



/([(/.],£)= / qm,^)d''r (All) 
J A 

is the force exerted by the fluid medium on the inclusion 
in the coarse-grained description. Thus, Eq. (AlO) gives 



(A12) 



which corresponds to the fundamental definition Eq. (2). 

Note that the force f{[(f)],£) is relevant at the coarse- 
graining level where the system is described by {[4>],f.). It 
is a force averaged over the microstates /i ^ ([(^],.^): in- 
deed, in the fundamental microstate fi, the force exerted 
by the medium on the inclusion is 



J A 



(A13) 



which verifies 

/ {q^)^^m,<^)d''r = m],i), (A14) 

J A 

where we have used Eq. (A2) and (All). 

If there are singularities in the energy density of the 
system at the boundaries of A, the mathematical def- 
inition f = Jj^q d'^r can become ambiguous. This is 
the case for instance when such singularities yield Dirac 
terms in q at the boundaries of A (see, e.g.. Appendix C). 
However, the force / must take into account all the terms 
that come from the presence of the inclusion, including 
boundary terms. Thus, in general, the integral over A 
has to be carried out using the following procedure: 



qd r = lim 



qd'^r, 



(A15) 



A, 



where A^ contains A and its boundary <S plus a shell of 
hypervolume e, so that each point of A and S is interior to 
Ac for all e > 0. This proccidurc; amounts to performing 
the integral infinitesimally outside the inclusion, in or- 
der to ensure that the whole inclusion is enclosed in our 
hypersurface of integration. Using this procedure does 
not change the result in the case where q is a piecewise 
continuous function. 

Note that we can also extract the hypervolume den- 
sity q of internal forces in the medium, at the coarse- 
graining level where the system is described by 
from Eq. (A9): 



(A16) 



Appendix B: Derivation of the stress tensor T of the 
fluid medium 



In this Appendix, we consider a fluid medium with 
Hamiltonian density 'H{(j),'^(j)) without any inclusion. 
The stress tensor associated to H can be derived from 
Noether's theorem (see, e.g., Ref. [30]). Indeed, it is the 
Noether current associated to the translation invariance 
of the system. It is obtained by considering the following 
infinitesimal translation: 



(r) 



r + Sr , 

(j){r) + d(j){r) 



(Bl) 



with constant Sr, and with (5^(r) = —'V4>{r) ■ 6r{r), i.e., 
S'^cp = 0. The stress tensor given directly by Noether's 
theorem is often called the canonical stress tensor (or 
more generally, in time-dependent field theories, the 
canonical energy-impulsion tensor). Indeed, it is its di- 
vergence that appears in the conservation law associated 
with translation invariance, so any tensor with vanishing 
divergence can be freely added to this canonical stress 
tensor. Such modifications of the stress tensor are often 
used in field theory, for instance to ensure its symme- 
try or its scale or conformal invariance [31]. However, 
in the case of perturbative embedded inclusions, where 
the stress tensor is defined everywhere in the system, 
even inside the inclusions (see Appendix C), Stokes' the- 
orem can be applied to the integral of the stress tensor 
on any closed hypersurface, so that the divergence-free 
"improvement" terms do not contribute to Casimir-like 
forces (this argument is used in Ref. [32]). We will not 
discuss these terms further. 

In the present work, our field-theoretic model describes 
a fluid medium where the interactions are supposed to 
be short-ranged. Thus, another way of constructing the 
stress tensor of the medium is to adapt the definition of 
continuum mechanics. Let us define the stress tensor T 
by the relation 



dfi=Tijnjd'^-'^r, 



(B2) 



where df is the infinitesimal force that one side of the 
medium (side 1) exerts onto the other side (side 2) 
through the hypersurface d'^~^r, and n denotes the nor- 
mal to this hypersurface directed toward side 1 [17, 33]. 
The stress tensor T of the medium can be determined 
from its energy H using the principle of virtual work. 
For this, we start by cutting (virtually) a subpart A of 
the medium, with energy 



Ha 



n(6,Vd))d''r 



(B3) 



Let us call B the rest of the medium (i.e., B = R'' \ ^), 
and S the interface between A and B. Let us now con- 
sider an infinitesimal transformation of the medium, as 
defined in Eq. (3). We consider that each fluid parti- 
cle keeps the same </> during the displacement: 6(j){r) = 
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— V(/)(r) • Sr{r), i.e., S'-^(f> = 0. The variation during this 
transformation of the energy of the cfosed system 
initially in A reads, at first order in e: 



-L 



d{dj(p 

dn 



did)5ri(fr 



dicp 



Tij 6ri d 



(B4) 



This variation of energy can be equated to the work 6W 
done during the transformation by the external forces 
acting on the closed system initially in A: 



SW= [ w,6rid^r+ [ Tijuj Snd^-'^r , 
J A JS 



(B5) 



where w is the hypersurface density of forces exerted by 
the exterior on the fluid medium, so that the integral 
on A represents the forces exerted by the exterior of the 
medium on the closed system initially in A. Meanwhile, 
the integral on <S represents the force exerted by the rest 
of the medium on the closed system initially in A, which 
has been expressed using the stress tensor thanks to its 
definition Eq. (B3). Since the energy balance SH = 6W 
must hold for any virtual deformation field Sr, we can 
identify the stress tensor of the medium: 



(B6) 



In addition, the identification of the bulk term gives the 
hypervolume density q of internal forces in the medium. 
Neglecting inertia, which is possible if the timescales con- 
sidered are sufficiently large, Newton's second law ap- 
plied to each fluid particle of the medium gives q = —w, 
where w represents the hypervolume density of external 
forces (see Appendix A). Thus, we obtain: 



-Wi 



dn 



dicj) = djTij. (B7) 



This derivation of the stress tensor shows that the 
stress tensor itself is built by assuming that each fluid 
particle keeps the same (j) during a displacement. It is 
therefore not surprising that we find the stress tensor 
when we calculate the force /^^^ under this assumption. 

With our mechanical definition of the stress tensor, 
we have obtained, from the principle of virtual work, a 
stress tensor T that is identical to the canonical stress 
tensor given by Noether's theorem. In field theory, the 
stress tensor T'^ generally written is in fact T[j = —Tji 
[30, 31], but this apparent difference is just a matter of 
convention. Note that our stress tensor T is fully defined 
(not up to a term with vanishing divergence) because we 
have asked it to give the force exchanged through any 
infinitesimal hypersurface. 



Appendix C: Stress tensor T' of the fluid medium 
with a perturbative embedded inclusion 

Let us now consider the fluid medium with a perturba- 
tive embedded inclusion in it: its effective Hamiltonian H 
corresponds to Eq. (17). Carrying out the same reason- 
ing as in the previous section, using a generic infinitesi- 
mal deformation 6r, enables to identify the stress tensor 
T' of the composite medium comprising the perturbative 
inclusion: 



T'. 
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{n + viA)Si, 



dn 

d{d,c^) 



dicj). 



(CI) 



The same reasoning also enables to identify the hy- 
pervolume density of internal forces at each point of the 
composite medium as q[ = djT-j. Explicitly, it gives 



SH 

— V<j) + VVlA, 



(C2) 



where H is defined in Eq. (17), and its functional deriva- 
tive with respect to (p is given by Eq. (19). Note that the 
gradient of 1^, and thus q' , features Dirac singularities 
on the contour of A, i.e., on <S. 



Appendix D: Pointlike embedded inclusion or 
non-embedded influencing object 

Let us consider a pointlike object in r = ^ that is 
coupled to ^{i): 

H = [ ni(i),V(p)d''r + v{4>{e)) 

= / [n + V6{r - £)] d'^r . (Dl) 

Here, the total variation of H during the infinitesimal 
transformation (3) reads, at first order: 



5H 



-J 

Jm 



dV 





( \] 







'^d\ 



{<!>{£)) [5(1) + V4> ■ 5i] . 



Thus, the first route, which corresponds to 
yields 

dV 



(D2) 
= 0, 

(D3) 



Meanwhile, if we follow the second route, i.e., if 6^(j) 
5(l) + V<f)-5r = 0, Eq. (D2) becomes 



SH = - 



f 
I 



dn_ f dn 

djTij Sn d'^r , 



d^(j) Sri d'^r 



(D4) 
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where wc have used Eq. (10). In spite of this simple 
expression, we must remember that the energy density 
H + V5{r — £) has a singularity in £. To deal with it, let 
us write 



6H = -lim [ djTijdndf-r 




where denotes the hyperball of radius e centered on 



I. This relation enables to identify /^^^ as 
= iim / djTijd'^rei 

= lim / TijHjdf'-'^rei, (D6) 

where iSf denotes the hypersphere of radius e centered on 
I and n is its exterior normal. 

Thus, the difference between the two ways of keeping cj) 
constant while varying I remains for pointlike objects. As 
in the case of extended objects, the force /^^^ is adapted 
to a pointlikc embedded inclusion, while the force f^^'^ is 
adapted to a non-embedded pointlike influencing object. 
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